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ABSTRACT 

We use the Bullet Cluster (1E0657-56) to investigate the extent to which star formation in cluster 
galaxies is influenced by ram pressure from supersonic gas (Mach 3) during a cluster merger. While 
the effects of ram pressure have been studied for individual galaxies infalling into galaxy clusters, this 
system provides a unique opportunity to investigate the impact of dramatic merger events on the 
cluster galaxy population. In this analysis we use Spitzer IRAC data to study star formation. At the 
rcdshift of the cluster the 6.2 (im PAH feature is redshifted into the 8 (im band, enabling use of the 
1114.5-iiig color as a proxy for specific star formation rate. We find that the color distribution on the 
two sides of the shock differ by less than 2a, and conclude that ram pressure from the shock front has 
no dramatic, immediate impact on the star formation of cluster galaxies in the Bullet Cluster. 

Subject headings: galaxies: clusters: individual (1E0657-56, The Bullet Cluster) - galaxies: evolution 



1. INTRODUCTION 

It has long been observed that there is a co rrelation 
betw een galaxy morphology and local density (|Dresslerl 
1980), known as the density-morphology relation. Al- 
though there are various physical mechanisms that can 
transform star-forming late-type galaxies into quiescent 
early-type galaxies, it is unknown whether this trans- 
formation occurs mostly in the galaxy group or galaxy 
cluster environment. 

Preprocessing of galaxie s in the group environment 
has been argued by some (jZabludoff fc MiilchaevllTQM 
iKodama et al.lr200lD to be a principal mechanism driv- 
ing the predominance of the passive galaxy populations 
observed in galaxy clusters. In galaxy groups, the veloc- 
ity dispersions are sufficiently low that processes such 
as galaxy-galaxy interactions can quench star forma- 
tion and morphologically a lter galaxies. In addi t ion to 
galaxy-galaxy interactions, iKawata fc Mulchaevl (2008) 
show that "strangulation" is an important mechanism 
that leads to quenched star formation in the group envi- 
ronm ent. Using a cosmological c hemodynamical simula- 
tion, IKawata fc Mulchaevl (|2008h show that most of the 
hot gas in a disk galaxy is stripped away, thus cutting 
off a source of new cold gas (some of the hot gas cools 
to form cold gas), which is necessary to maintain star 
formation. 

Others have argued that cluster-specific processes are 
responsible for the high e arly-t ype fraction. The sim- 
ulations of iBerrier et al.l ((2008) argue that the trans- 
formation of late-type to early-type galaxies is mostly 
attributed to processes internal to the galaxy clus- 
ter environment, rather than the group environment. 
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IBerrier et al.l (|2008h find that most of the galaxies in a 
galaxy cluster are accreted directly from the field, rather 
than from infalling galaxy groups. In this case, mech- 
anisms such as galaxy har assment dMoore et all |1996|) 
and ram pressure stripping (jGunn fc Gottlll972[) that are 
most efficient in massive clusters must drive this trans- 
formation. 

Simulations have shown that ram pressure stripping 
can remove some or all of the interstellar medium 
(ISM) fro m late-type galaxies as the y enter the c luster 
potential dRoediger fc Henslerl [200l iMaver et~ail 120061: 
Briiggen fc De Luciall2008D . Studies of galaxies in clus- 
ters with observed HI deficiencies or truncated disks sup- 
port the idea that ram pressure stripping can transform 
a ga s-rich spiral galaxy into a n anemic spiral galaxy 
(e.g. IVollmer et all 120011 . 120081) . which may ev entually 
evolve into a lent i cular (e.g. iBekki et al.l [20021 . but see 
iBoselli fc Gavazzil ()2006D ). Ram pressure has also been 
sho wn to trigger bursts o f star f ormation in cluster galax- 
ies dGavazzi et aL I2001L |2003| ). and recent simulations 
(|Kronberger et a" I2008D predict enhanced star formation 
rates in galaxies undergoing a ram pressure event. 

While previous studies have examined the impact of 
ram pressure on individual galax i es falling into the clus- 
ter potential (jCortese et al.ll2007t ICtiung et al.ll2007f ) the 
Bullet Cluster provides a unique opportunity to examine 
the effect of ram pressure induced by cluster mergers. 
Although galaxies in any cluster environment undergo 
ram pressure as they travel through the ICM, galaxies 
in a major merger event such as those we observe in the 
Bullet Cluster will experience a dramatic enhancement 
in the pressure since P tx V 2 . It is thus possible that 
this brief transient phase may have a significant impact 
on the properties of the cluster galaxy population. The 
Bullet Cluster is an ideal site in which to quantify the 
importance of such merger-induced ram pressure. 

In this paper we conduct an initial exploration of the 
impact of the shock front upon polycyclic aromatic hy- 
drocarbon (PAH) emission, and hence star formation, 
in cluster galaxies. We analyze Spitzer data taken 
with the Infrared Array Camera (IRAC; iFazio et al.l 
l2004f) . in conjunction with optical spectra observed with 
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the Inam ori Magellan Area! Camera and Spectrograph 
(IMACSi lBieelow et al.lll998h . 

2. GEOMETRY AND CLUSTER PROPERTIES 

The Bullet Cluster, at a redshift of z = 0.296, is 
a head-on collision between two clusters of unequal 
mass, which we c all th e main cluster and subcluster. 
iMarkevitch et alj (|2002D [M02] observed a sharply de- 
fined bow shock front, driven by the subcluster and prop- 
agating in the gas of the main cluster with a velocity of 
4740±™ km s" 1 (Mach number M = 3.0±g;||) (Marke- 
vitch et al. 2008, in preparation). Hydrodynamic simula- 
tions showed that the subcluster itself has a lower veloc- 
ity relative to the main cluster, ~ 2700 km s _1 (Springel 
& Farrar 2007; see also Milosavljevic et al. 2007). While 
the subcluster and the shock front might be expected to 
move together, the gravity of the subcluster causes the 
main cluster's gas in front of the shock to flow onto the 
subcluster. As a result, the shock front has a higher ve- 
locity in the reference frame o f that gas inflow, which 
is what is measured in X-rays. IMarkevitch et all (|2004l ) 
constrain the inclination angle of the collision, with re- 
spect to the plane of the sky, to i < 8°, based on the Mach 
number of the shock and subcluster relative l ine-of-sight 
velocity (-600 km s" 1 ) (jBarrena et alJl2002f >. From the 
velocity and geometry of the system, it is computed that 
the subcluster exited the main cluster core ~0.15 Gyr 
ago. 

The simple geometry of the Bullet Cluster makes it 
an ideal system in which to study how galaxy properties 
such as star formation are affected by ram pressure. Be- 
cause the merger is nearly in the plane of the sky, the ve- 
locity of galaxies with respect to the intracluster medium 
(ICM) is well constrained, providing a good quantitative 
measure of the ram pressure exerted by the gas on the 
galaxies. In addition, the recentness of the collision (0.15 
Gyr) allows us to use sensitive indicators of recent star 
formation, direct and indirect measures of the ionizing 
flux from O/B stars with lifetimes of 0.1 Gyr. 

3. PHYSICAL ASSUMPTIONS 

The analysis in this paper requires several assumptions 
related to the star formation rate and the three-space 
galaxy distributions for the main cluster and subcluster. 

3.1. Star Formation Rate 

To probe star formation rates, we use the PAH emis- 
sion at rest wavelength A = 6.2/im, which at the clus- 
ter redshift is observed in the 8/xm IRAC band. Ultra- 
violet (UV) photons that are produced by young hot stars 
are known to excite PAH molecules whi ch then emit in 
the A = 6.2/nm and other infra red bands (jLeger fc Pugetl 
ll984tlAl"iamandola et al.lll985l ). More recent studies have 
shown that the 6.2/mi luminosity is correlated with to- 
tal in frared luminosity (|Peeters et al. 2004; Brand l et al. 
2006) , which is a tracer of star formation rate (jKennicutt 
1998). However, PAH emission is not a unique signa- 
ture of star fo rmation, which can also be excited by vis- 
ible photons ()Uchida et all 119981 : iLi fc Drainell2002D. as 
well a s by active galactic nuclei (AGN) |Freudling et alJ 
12003ft . It has also been sh own that PAH emission 
is dependent on metallicity (|Engelbracht et all 120051 : 
Smit h~et al.ll2007T) and thus not a robust star formation 



rate indicator for galaxy samples that span a large range 
in metallicity (jBoselli et al.ll2004h . 

For the Bullet Cluster, both effects are expected to 
be second order. Our data are negligibly affected by 
AGN contamination (see §5), and our sample is domi- 
nated by massive galaxies, for which the metallicity vari- 
ations are not dramatic. Moreover, the metallicity dis- 
tribution should be the same for galaxies on both sides 
of the shock. The PAH flux should nevertheless be con- 
sidered a rough proxy for star formation rate rather than 
a precision tracer. 

To normalize the 8/iin flux, we use the 4.5/im flux as a 
proxy for mass. The 4.5/im flux originates mostly from 
the old red giant stars, and probes the stellar mass of a 
galaxy at the cluster redshift. By using the ratio of 8/im 
to 4.5/xm flux or 1114.5-1118 color, we can roughly trace 
star formation rate per stellar mass. Specifically, 8/im 
emission in excess of that expected from the Rayleigh- 
Jeans tail of the cold stellar component is assumed to 
arise from star- format ion induced PAH emission. We 
estimate the level of impact that the M~3 shock has 
had on star formation activity in individual galaxies by 
looking at the 1114 5-mg color as a function of distance 
from the shock front. 

3.2. Galaxy Distribution and Projection 

We proceed with our analysis under the simplifying as- 
sumption that galaxies that lie ahead of the shock front 
are main cluster "pre-shock" galaxies, and those behind 
the shock front are main cluster "post-shock" galaxies 
(Figure [T]). The post-shock galaxies have already expe- 
rienced the ram pressure exerted by the gas behind the 
shock, whereas the pre-shock galaxies have not yet been 
affected. The true physical situation is clearly more com- 
plex, with contamination from the subcluster galaxies 
as well as background/foreground sources not associated 
with the Bullet Cluster. Although there is some contam- 
ination from the subcluster galaxies, we assume that the 
galaxy populations on both sides of the shock front are 
dominated by main cluster galaxies based on estimates of 
the total stellar mass. The stellar mass of the main clus- 
ter is approximately an order of magnitude larger than 
the stellar mass of the subcluster. 

We also note that any effect of contamination from the 
subcluster galaxies will be strongest at the shock front, 
which bisects the subcluster brightest group galaxy. Sub- 
cluster galaxies just ahead of the shock front experience a 
strong ram pressure as they collide with the main cluster 
ICM at ~4700 km s^ 1 , whereas subcluster galaxies just 
behind the shock are somewhat shielded by the shock 
front and feel a much weaker ram pressure with a slower 
gas inflow velocity of ~1600 km s _1 . Therefore any effect 
of ram pressure on star formation rate has the opposite 
sign for subcluster galaxies near the shock front, com- 
pared to the main cluster galaxies. 

We assume that there is a sharp division between post- 
shock and pre-shock galaxies, where all galaxies "behind" 
the shock front (leftward of the curve in Figure [lj have 
experienced the effects of the shock. However, because of 
the projection of galaxies and curvature of the shock front 
surface along the line of sight, there is a gradual transi- 
tion from post-shock to pre-shock galaxies. This projec- 
tion effect softens any step-like increase/decrease in color 
(or star formation rate) across the shock boundary. By 
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generating galaxies in a random distributi on according to 
the N avarro, Frenk, and White (NFW) (jNavarro et al.1 
|1997|) profile, and assuming the shock front is an ax- 
isymmetric cone projected in the plane of the sky and 
spanning —550 kpc from the axis of symmetry, we es- 
timate that the transition region during which the ob- 
served galaxies go from being fully post-shock to fully 
pre-shock galaxies is ~0.5 Mpc. 

4. OBSERVATIONS 

4.1. IRAC 

Observations of 1E0657-56 were taken on December 17, 
2004, with Spitzer/IRAC. Data were collected from all 
four IRAC channels - 3.6//m, 4.5//m, 5.8/im, and 8^tm, 
in full array readout mode. A medium scale, cycling 
dither pattern was used, with 20 pointings. At 2 frames 
per pointing and 100 s exposure time per frame, the to- 
tal integration time per IRAC band was 4000 s. The 
native IRAC pixel scale is 1.22 arcsec, but the final re- 
duced images were set to have a pixel scale of 0.86 arcsec 
pixel -1 , with an effective field of view (FOV) of 3.7 x 3.7 
arcmin covered by the four channels. The frames were 
processed using the Spitzer Science Center (SSC) IRAC 
Pipeline, and mosaics created from the basic calibrated 
data (BCD) frames using MOPEX. 

4.2. IMACS 

We analyzed optical spectroscopic data taken with 
Magellan/IMACS in 2005 and 2006. In 2005, 430 ob- 
jects were observed with four masks, with exposure times 
ranging from 4.8 to 5.4 ks. During this run, red sequence 
galaxies were given highest priority as targets. The fol- 
lowing year, two additional masks were obtained, with 
exposure times of 10.8 and 12.0 ks. Target priority was 
given to objects with relatively red 1114.5-mg colors. The 
wavelength range of all spectra are 4000 A— 9000 A, with 
a dispersion of 2 A per pixel. We note that spectroscopic 
selection criteria for galaxies on both sides of the shock 
front are identical. 

5. DATA ANALYSIS 

5.1. IRAC Photometry 

Source detection and photometry of the IRAC images 
are performed with Source Extractor (Bertin & Arnouts 
1996). Sources are identified in the 4.5/um image with 
a 3.5a detection threshold. Magnitudes are measured 
within a 5" aperture. We choose to apply no aperture 
corrections to the photometry. In this analysis we are 
primarily interested in colors rather than total magni- 
tudes, which for mg-ms.g and ni4. 5-1x13. 6 would change 
by only —0.05 magnitudes with the corrections applied. 
Since the corrections themselves are uncertain by up to 
10% ( [http://ssc. spitzer.caltech.edu/irac/calib/extcal), 
excluding systematics, we use raw aperture magnitudes 
for all analyses. In addition, aperture corrections are 
expected to be minor when examining colors rather than 
absolute magnitudes. 

5.2. IRAC Sample Selection 

The curvature of the shock contour indicates that the 
direction of propagation lies within —5 degrees of the 
east-west axis (M02). The shock front was defined as 



the X-ray brightness contour whose most westward point 
lies at 104.6 degrees in right ascension (Figure [l}, with 
a maximum and minimum declination of -55.908 and - 
55.985 degrees. Objects located above or below (north or 
south) the shock front were not included in our analyses, 
since these sources are not in the direct path of the shock 
front and thus should be minimally affected by the prop- 
agating shock. We limit our analyses to the region of 
spatial overlap between the four IRAC bands, detecting 
a total of 758 sources within this "analysis region" . 

We eliminate unsaturated stars fro m the IRAC catalo g 
using Magellan optical photometry (| Clowe 
In the inner region of the R-band image where the point 
spread function (PSF) is small (-0.45 arcsec FWHM), 
objects with a half-light radius of 2.6 to 3.2 pixels (0.111 
arcsec/pixel) and an R-band aperture magnitude of 18 to 
24.1 (with a 7.8 pixel aperture), are flagged as stars. Due 
to a gradient of the PSF in the image, half-light radius 
and aperture magnitude cuts are scaled accordingly to 
provide a catalog of unsaturated stars. Saturated stars 
are identified manually by inspecting the Magellan R- 
band image to look for saturated objects. 

To minimize contamination from bright foreground 
galaxies, we impose a magnitude limit of 1114.5 > 14.20 
- the magnitude of the brightest cluster galaxy (BCG) 
in the main cluster. We also apply a magnitude cut of 
1114.5 < 17, which is sufficiently faint to detect most of the 
8/xm sources. While going deeper to m 4 . 5 < 18 increases 
the number of cluster candidates, we opted for a more 
conservative 1114.5 < 17 magnitude limit so as to mini- 
mize the background contamination. Most importantly, 
the analyses described in this paper produce the same 
result within error bars, whether we apply a magnitude 
limit of 1114.5 < 17 or 1114.5 < 18. 

Of the 758 sources, 18 are brighter than the main clus- 
ter BCG, 176 are stars, leaving us with 564 sources. 
Eliminating sources fainter than 1114.5 — 17.0, we are left 
with 169 sources in the direct path of the shock front 
and 31 that are above/below the shock front. The dis- 
tance from likely star-forming galaxies to the shock front 
is illustrated by arrows in Figure [1] 

5.3. IMACS Spectroscopy 

To reduce contamination from background and fore- 
ground sources, we reduce a total of 458 IMACS spec- 
tra with IRAF and IDL, then cross-correlate them with 
template spectra, using the IRAF task xcsa o to obtain 
redshi fts. We use the same templates as iTran et all 
(2005), each spectrum representing a different type of 
galaxy - giant elliptical, E+A, Sb spiral, and an emis- 
sion line galaxy (ELG). Of the 458 spectra, we obtain 
reliable redshifts for 326 sources. W e augment our spec- 
troscopic data set with 133 redshifts (|Barrena et al.ll2002l . 
priv. comm), which more densely sample the cluster 
core, yielding a total of 459 sources with known redshifts. 

We next associate the spectroscopic redshifts with 
IRAC sources for the 114 spectra that lie within the 
IRAC field-of-view. Of these 114 sources, 67 are cluster 
members with redshifts that lie within ±2000 km s _1 of 
the mean cluster redshift (z — 0.296), and 47 are inter- 
lopers. Sixty-three out of the 67 cluster members, and 
42 out of the 47 interlopers, meet the 1114.5 < 17 criteria. 
This leaves us with 134 "cluster candidates", which in- 
clude all IRAC sources that meet our magnitude criteria, 
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are located within the direct path of the shock front, and 
exclude all known interlopers. 

6. RESULTS AND DISCUSSION 

6.1. Impact of the Shock 

We use the ni4.5-ni8 color to probe relative star forma- 
tion rates per unit stellar mass. At the redshift of the 
cluster (z = 0.296), the prominent PAH emission band at 
rest wavelength 6.2/^m falls into the 8.0/im IRAC band. 
PAH emission, which is excited primarily by UV photons, 
most often originates fro m photodissqciatiq n regions of 
star formation sites (e.g. iHowell et aLll2007f l. Although 
PAH emission is som etimes also associated with AGN 
(|Freudhng et al.ll2003l ). only 2 out of 134 cluster candi- 
dates in the direct path of t he shock front have IRAC col- 
ors w ithin the AGN wedge (jLacv et al.ll2004t [Stern et all 
2005). Both AGN candidates have m4. 5 -m 8 < 0.6, and 
thus do not contribute spurious star formation signa- 
tures. 

Figure [5] shows color as a function of distance from the 
shock front for 134 cluster candidates. The dotted hori- 
zontal lines show the exp ected color of an e lliptical and 
Sbc type galaxy at z=0.3 (jAssef et al.ll2008l ). These pre- 
dicted colors come fro m a set of low-reso lution spectral 
templates, derived by lAssef et al.l |2P08), using optical 
and near-infrared photometry of over 16,000 galaxies in 
the NOAO Deep Wide-Field Survey Bootes region, and 
redshifts from the AGN and Galaxy Evolution Survey 
(AGES). For elliptical galaxies the 8/im emission comes 
from the Rayleigh- Jeans tail of the ste llar component, fo r 

which the m 4 . 5 -m 8 0.1 at z=0.3 (jAssef et al.ll2008[ ). 

Colors redder than the elliptical locus are indicative of 
PAH emission in the 8/im band. 

A majority of the IRAC sources have colors consistent 
with red sequence galaxies in the cluster - only 17% of 
the cluster candidates have a color ni4.5-m8 > 0.5, and 
6% of the sources have m 4 . 5 -m 8 > 1. Although there is 
a selection effect in favor of red sequence galaxies with 
the spectroscopic data, Figure [2] uses the entire IRAC 
sample, for which there is no bias towards quiescent el- 
lipticals. 

The spectroscopic data is useful for eliminating in- 
terlopers, but the sample of confirmed cluster members 
(solid points in Figure[2]& Figure[3j) is of limited value in 
the current analysis because the spectroscopic program 
preferentially targeted red sequence galaxies. Indeed, 
only 2 of the 23 cluster candidates with nx^-mg > 0.5, 
presently have spectra confirming their membership. 

Instead, we focus in the current analysis on the 
full sample of 134 candidate members and assess 
whether there is evidence for a correlation between the 
luminosity-weighted color of cluster members and their 
distance from the shock front. For each galaxy we com- 
pute the distance from the shock (Figure Q]) , and then 
compute the integrated color of all galaxies in 275 kpc 
wide bins (Figure [3]). For each bin we use the boot- 
strap technique with 10,000 realizations to calculate the 
uncertainties. The dotted lines correspond to the total 
integrated color of all objects on each side of the shock 
front, with the width of the shaded grey region corre- 
sponding to the la confidence interval calculated via the 
bootstrap method. 

For the individual data points on Figure [3l the size of 



the bootstrap error bars indicates that the mean color 
is strongly sensitive to the small subset of strongly star- 
forming galaxies. For example, in the bin directly ahead 
of the shock front, there are 2 (out of 22) galaxies that are 
responsible for the relatively high integrated color. One 
of these two galaxies is a spectroscopically confirmed, 
face-on spiral; the other is a large disk galaxy whose op- 
tical color is bluer than the red sequence. 

The horizontal dotted-dashed line in Figure [3] shows 
the predicted colo r of an elliptical galaxy at z = 0.3 
(jAssef et al.l I2008D . The data shown in Figure [3] are 
meant to show an excess of color in the post-shock and 
pre-shock galaxies, in comparison to what one would 
expect from a purely quiescent elliptical galaxy. Over- 
all, Figure [3] shows that there is no drastic change in 
color across the shock boundary. The integrated color of 
the post-shock versus pre-shock galaxies is consistent to 
within 2cr. 

The above technique is a rather blunt means of as- 
sessing the impact of the shock front, given that it in- 
cludes the red sequence galaxies, where we expect mini- 
mal change, as well as the star-forming galaxies. An al- 
ternate approach is to use the Kolmogorov-Smirnoff (KS) 
test to understand how statistically different the popula- 
tions are on either side of the shock. We first compare the 
full distributions on the two sides of the shock front. Us- 
ing all 134 candidate members, a KS-test indicates that 
the colors on the two sides of the shock differ at only 
the 1.3ct level; consequently, the data do not exclude the 
hypothesis that the galaxies are all drawn from the same 
population. We note that the full sample includes ellipti- 
cal galaxies, whose 1114 5-mg color simply reflects the cold 
stellar component of the Rayleigh- Jeans tail rather than 
star formation. All galaxies are included in our initial 
KS-test because we are looking for changes in the global 
color distribution of galaxy colors across the shock front. 

To avoid the gas-poor elliptical galaxies whose 104.5- 
mg color is not a good proxy for specific star formation 
rate, we can exclusively use galaxies whose n^.s-mg color 
is redder than some value. Comparing the color dis- 
tribution across the shock front for galaxies with 1114.5- 
mg > 0.2, m 4 .5-ni8 > 0.5, and 1114.5-nig > 1, the KS-test 
indicates that the colors on the two sides of the shock 
front differ at the 0.8(7, 1.5cr, and 1.2cr level, respectively. 
Although we avoid the elliptical galaxies by imposing 
color cuts, we also introduce an uncertainty since objects 
may move in and out of the sample due to the impact 
of the shock front, without leaving a known trace on the 
final color distribution. However, we emphasize that the 
KS-test for all of our color cuts shows no evidence for sig- 
nificant change in the color distribution from pre-shock 
to post-shock galaxies. 

For comparison, we also investigate the optical prop- 
erties of our IRAC sample. Figure [4] presents the color- 
magnitude relation (CMR) for the 88 objects out of the 
134 IRAC cluster candidates with R < 20.5, using our 
Magellan data and again excluding spectroscopically con- 
firmed interlopers. As indicated by the IRAC colors, 
Figure H] confirms that most of the IRAC sources are red 
sequence galaxies. However, to determine exactly which 
objects can be classified as part of the red sequence, we 
need three parameters of the CMR - the slope, zero- 
point, and dispersion. 

We adopt a CMR slope of -0.076, derived from 57 X- 
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ray clusters examined by lLopez-Cruz et al.l {2004) , with 
rcdshifts ranging from 0.02< z <0.18. To obtain the zero- 
point of the CMR relation, we fit the data with this fixed 
slope. This fit yields a zero-point of 3.927, which is ~0.1 
mag brighter than the expected zer o-point calculated 
from t he zero-point-redshift relation of lLopez-Cruz et al.1 
(|2004f ). Finally, we calculate the dispersion by fitting a 
Gaussian function to the CMR residuals, excluding out- 
liers. We derive a = 0.136 mag, indicating that our dis- 
persion is indeed dominated by photometric uncertainty. 

We classify all objects within 2a of the CMR as red 
sequence galaxies, indicated by the dotted lines in Fig- 
ure |U Out of the 7 objects with m 4 . 5 -m 8 > 0.5 and 
R < 20.5, 5 of them lie more than 2a below the CMR 
(with one object just barely within 2a). These red IRAC 
sources consitute 5/12 or 42% of the objects below the 
2a line. The optical data thus confirm that objects with 
red 1114.5-nig colors correspond to optically blue and likely 
star-forming galaxies. 

There are a few caveats in interpreting the results dis- 
cussed above. Foremost, there remains contamination 
by foreground and background sources for the galaxies 
that lack spectroscopic confirmation. Foreground inter- 
lopers among the star-forming galaxies may depress the 
observed significance of any color change across the shock 
front. It is unlikely however that removing these fore- 
grounds will create a significant difference in the mean 
m4.5-mg color across the shock boundary. A total of 55% 
of IRAC sources with color m^-nig > 0.5 are known in- 
terlopers, while the fraction of known interlopers is only 
22% for sources with 1114.5-1118 < 0.5. Because the frac- 
tion of interlopers is higher for late-type galaxies than 
early-type galaxies, eliminating these objects will likely 
cause the mean color on both sides of the shock front 
(Figure [3]) to approach the quiescent value. 

6.2. Physical Interpretation 

The IRAC data show no significant change in color 
across the shock boundary (Figure [3]), indicating that 
the gas behind the shock front does not have a large 
impact on the current or recent star formation activ- 
ity in the individual galaxies. One possible explanation 
is that ram pressure stripping may deplete existing gas 
reservoirs in the outskirts of galaxies without disturbing 
current/recent star formation sites that may be more 
centra lly located in the disk. Following iGunn 
(1972), a galaxy moving face-on through an ICM expe- 
riences a ram pressure of P = pV 2 . If ram pressure 
exceeds the galaxy's internal gravitational pressure, ram 
pressure stripping occurs. The gravitational pressure can 

A/ 2 

be expressed as P = 2-kGx — sj f 2L , where G is the gravi- 

tational constant, x is the interstellar gas to stellar mass 
ratio, M s tar is the stellar mass, and r g is the galaxy ra- 
dius. A galaxy with stellar mass M* = 10 10 Mq and 
x = 0.6, should experience ram pressure stripping at 
a distance of ~4 kpc from its center, assuming that it 
moves face-on through an ICM of density 4 x 10~ 3 cm" 3 , 
at a velocity of ~2700 km s _1 . This is somewhat beyond 
the expected disk sc ale length for a gala xy with stellar 
mass M* = 10 10 M Q (jPizagno et al.ll2005t ). Therefore, we 
might not expect to see evidence of disturbed star for- 
mation due to ram pressure, since much of the current 
star formation occurs within the disk scale length. How- 



ever, the gas density and galaxy velocity (with respect to 
the ICM) are highly uncertain at any significant distance 
from the shock front, since both quanities change with 
time and position in the cluster. 

If ram pressure stripping destroys gas reservoirs in 
galaxies that host significant star formation, we might 
expect to see an excess of post-starburst galaxies in ~10 
Myr (after a galaxy crosses the shock front), once the 
O stars no longer significantly contribute to the stellar 
population. A timescale of 10 Myr corresponds roughly 
to 75 kpc in Figure [3l On such short timescales, it 
is difficult to detect an excess of any particular galaxy 
type because of small number statistics (i.e. -there are 
not enough galaxies within a single 75 kpc bin). Also, 
post-starburst galaxies are difficult to distinguish from 
quiescent elliptical galaxies with phot ometry, since bot h 
populations lie near the optical CMR (jTran et al J 1 20071) . 

However, t he lifetime of a post -starburst galaxy is 
~1 Gyr (e.g. Qu intero et al.l I2004T ). providing us with 
a longer timescale (or distance from the shock front) 
in which we can search for an excess of these galaxy 
types. In future work, we will spectroscopically iden- 
tify the E+A galaxies in our sample. A higher fraction 
of E+A galaxies behind the shock front versus ahead of 
the shock front, coupled with no change in current star 
formation rates across the shock boundary, would favor 
ram pressure stripping of gas reservoirs, rather than some 
effect on current star formation. 

Although Figure [3] shows that there is no significant 
change in color in the pre-shock versus post-shock galax- 
ies, the geometric effects and contamination from non- 
cluster sources may dilute the signal. Future work in- 
cludes analysis of IRAC and MIPS 24^m data observed 
over a wider area, which coupled with spectroscopic red- 
shifts will be less sensitive to geometric effects and con- 
tamination. However, even with our current level of con- 
tamination, it is clear that the ram pressure exerted from 
the supersonic gas does not have a dramatic impact on 
current star formation activity in the main cluster galax- 
ies. 

Another source of signal dilution is the large frac- 
tion of early-type galaxies in our sample, for which 
there are two possible explanations. First, galaxy pre- 
processing in the group environment may have already 
done its work on both the cluster and subcluster galax- 
ies, explaining the observed lack of late-type g alaxies 
HZabludoff fc Mulchaevl[T998t ICortese et al.l 120061) . Sec- 
ond, the galaxies we observe may have been affected 
by the cluster environment in such a manner that they 
have already convert ed th eir gas reservoi r into stars. 
iMarcillac etail (|2007l) and iBai et alj (|2007t) have inves- 
tigated star formation of galaxies in merging clusters at 
z ^0.83. They find evidence of triggered star forma- 
tion in inf ailing galaxies and galaxy groups, which cou- 
pled with mechanisms such as ram pressure stripping or 
galaxy harassment, can act to eventually quench star for- 
mation. In this scenario, many of the initially gas-rich, 
star forming galaxies in the main cluster would have con- 
sumed much of their available gas in an earlier epoch of 
triggered star formation, explaining the lack of late-type 
galaxies in our current sample. However, the relative lack 
of late-type galaxies cannot explain why we do not see a 
signal across the shock front, especially considering that 
even among the late-type population, we do not detect 
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a significant difference in color in the post-shock versus 
pre-shock galaxies. 

7. SUMMARY 

To trace specific star formation rate, we use IRAC color 
m4.5-ni8. The 8fim flux traces star formation via the 
PAH emission feature at rest wavelength 6.2/xm, which 
is redshifted into the 8/zm band at the cluster redshift. 
The 4.5^.m flux can be taken as a proxy for stellar mass. 

Using spectroscopic redshifts, we confirmed 63 out of 
200 IRAC sources (169 in direct path of shock front + 
31 above/below the shock front) as cluster members, and 
found 42 IRAC sources that are interlopers. Applying a 
magnitude cut at 1114.5 = 17, we are left with 158 cluster 
candidates, of which 134 are in the direct path of the 
shock front. 

We observe no significant trend in m4.5-mg as a func- 
tion of distance from the shock front. The KS-test re- 
veals that the sample of galaxies behind the shock ver- 
sus ahead of the shock front only differs in color by less 
than 1.5(7, for galaxies with various cuts in m4.5-ms color. 
This indicates that even in a dramatic merger event like 
the Bullet Cluster, the ram pressure induced from a su- 
personic collision does not drastically trigger or quench 
current /recent star formation. A possible explanation is 
that the ram pressure can remove the gas reservoirs in the 
outer disk, while preserving star formation in the central 
disk of a galaxy. While geometric effects, contamination 
from interlopers, and the lack of late- type galaxies in the 
Bullet Cluster could contribute to diluting a signal, they 
are unlikely to explain the complete lack of signal we ob- 
serve. Even among our small sample of late- type galax- 
ies, we detect no change in color between the pre-shock 
and post-shock galaxies to within 1.5a. 
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Fig. 1. — Bullet Cluster candidates (filled circles) shown with the shock front (a corresponding contour of constant X-ray surface 
brightness). Arrows indicate the distance between the object and the shock front for sources with IRAC color 1114 5-mg > 0.5. 
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Fig. 2. — The nu.g-mg color as a function of distance from the shock front for 134 cluster candidates. The lower and upper axes indicate 
projected distance in degrees and Mpc, respectively. Spectroscopically confirmed cluster members are filled in boxes. The vertical dotted 
line shows the location of the shock front. Sources to the left of this line are considered post-shock galaxies, and those to the right are 
considered pre-shock galaxies. The horizontal dotted-das hed lines a t m4.5-ms=-0.13 and m4.5-ms=2.43 represent the expected color of an 
Elliptical and Sbc galaxy at z=0.3, respectively (Assef et al. 2008). The error bar in the upper left corner shows the mean color error of 
the 134 candidates. 




Fig. 3. — The 1114. 5-ms color as a function of distance from shock front for cluster candidates in bins of 275 kpc, with the number of 
sources per bin denoted. Individual error bars are derived from bootstrapping method. Grey shaded regions show the integrated color for 
all objects behind and ahead of the shock front, ±1<7 derived from bootstrapping. 
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Fig. 4. — The optical color-magnitude diagram shown for 88 IRAC sele c ted cl uster candidates. The solid line is the color-magnitude 
relation whose slope and zero-point were adopted from Lopcz-Cruz et al. (2004). The dotted lines signify 2a from the fit, where the 
dispersion is obtained by fitting a Gaussian to the residuals. Objects overplotted with a cross symbol have IRAC color m4.5-m8 > 0.5 and 
are thus likely to be star-forming galaxies. 



